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7'he Junction of the angle-error detector is to provide 'pointing-error 
signals to the ground antenna control system, which allows operation in 
the autotrack mode once the satellite beacon has been acquired. The limita- 
tions on the accuracy of this system imposed by noise, phase jitter and 
Doppler ejjecis are evaluated and the optimum, design in terms of minimum 
mean-square error is developed. Design examples are given for both the horn- 
reflector antenna autotrack system and the precision tracker antenna system. 

I. INTRODUCTION AND SUMMARY 

To insure the at^iuisition and accnratc tracking of the Telstar com- 
mimication satellites, a .sequence of tracking modes is provided at 
the ground stations hi Anilo\'er, Maine and Pleumeur-Bodou, l-'rancc.^ 

Initial pointing directions to both the precision tracker^ and the 
horn-reflector antennas are provided from orbital data appropriately 
processed and up-dated for each satellite pass. Once the precision tracker 
acquires and tracks the satellite, the horn-reflector antenna can use the 
pointing directions received from the precision tracker control system 
to acquire the satellite beacon signal m its narrow l)eam\vidth. Finally, 
the autotrack system^ provides closed-loop automatic control of the 
horn-reflector antenna using eiior signals derived from the satellite 
beacon.* 



* After orbital data becomes Bufficiently accurate, it is possible for the horn 
antenna to acquire the satellite from initial pointing directions and then go 
directly into the autotrack mode without using the precision tracker. 
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The detection of these error signals is accomplished in the system 
described in this paper. The inputs to this system are obtained by 
means of a mode separation technique' in the waveguide of the horn 
antenna, one mode having a peak amplitude on target, the other having 
a null, similar to the sum and difference signals in conventional mono- 
pulse tracking systems. The characteristics of these input signals and 
noise are discussed in Section II. An analysis of the phase-lock de- 
tection scheme which converts these inputs into the desired antenna 
pointing-error signals is given in Section III. The accuracy of these 
pointing-error signals is shown to be critically dependent upon the 
degree of phase coherence achieved by the phase-lock loop, which is 
discussed in Section IV. The design of the phase-lock loop to minimize 
the mean-square phase error in the output signals is considered in 
Section V. 

In Section VI a numerical example is given for the optimum design 
of the phase-lock detector in the vernier autotrack system. Since the 
precision tracker also uses essentially the same angle-error detection 
scheme, a parallel design of this system is included for comparison. 

II. INPUT SIGNAL AND NOISE CHARACTERISTICS 

The function of the angle error detector is to develop electrical error 
.signals proportional to the pointing angle error, |S, between the antenna 
boresight and the actual satellite position. The expressions for the 
desired output error signals are 

tx = 8 cos u> 

Cy = P Sin if) 

where x and y are Cartesian coordinates in the plane normal to the 
antenna boresight (electrical) and ip is the angle which the projection 
on the x-y plane of the radius vector, R, to the satellite makes with the 
X-axis, (see Fig. 1). 

The information on the parameters ^ and <fi necessary for the error 
signals (1) is contained in the amplitude and phase, respectively, of 
the difference channel received signal relative to the sum channel* 
received signal. For a pointing error, /3, which is within the beamwidth 



•The designations "difference" and "sum" channels are carried over from 
conventional monopulse usage. For the autotrack system, these terms should be 
"TMio" and "TEu" channels, respectively. This analysis is applicable to a single 
plane in conventional (linearly polarized) monopulse if the angle ip is taken to be 
or 180 degrees. 
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Fig. 1 — Pointing error, j9, and the angle ip in the x-y plane which determine 
the orthogonal error signals: ti = & cos ip; e^ = ^ sin p. 

of the sum pattern (the TEn mode pattern in the horn-reflector an- 
tenna), the received signals in the sum and difference channels can be 
expressed as^ 

eM) = E,(R,^) cos U,t + 0.(0) + N,{i) 

(2) 
ed(n = ■n^E.iRS) cos (wfi/ + did) + ^) + Ndit) 



where 



6)b = 2ir X frequency of satellite beacon transmitter, 
R = range of the satellite, 
E^{R,fi) = sum channel signal amplitude, 

Tj = difference channel relative sensitivity, 



1 

El 



HE, 



A/3 



^=u 



and 



Oi = e,(l) + djt), is the signal phase relative to a reference phase, 
dr = 0, plus a random phase fluctuation, 6„(t), discussed 
below. 
N,{t) and Nd(t) are the thermal noise components at the inputs of 
the sum and difference channels, respectively, whose one-sided power- 
spectral densities are assumed identical and equal to 



where 



Oat = kTeg watts/cps 



k = 1.38 X 10"'' watt-sec/°K 
Teg = equivalent receiver noise temperature, °K. 



(3) 
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The random phase fluctuation, 5„(0, results from the frequency 
instability of the various oscillators in the system, principally the 
beacon oscillator in the satellite, since elaborate frequency stabilizing 
techniques are not feasible from weight and space considerations. The 
one-sided spectral density of the resultant phase fluctuation can be 
expressed as 

$e = - — ^ radVcps (4) 

where 

Tee = equivalent coherence time of the system oscillators.* 
For the purpose of this analysis, the thermal noise terms in (2) will 
be represented by the usual in-phase and quadrature notation*'* 

N{t) = X(t) cos (oii,t + di) + F(n sin (utt + tf.) 

where Xit), Yit) = independent Gaussian random voltages with one- 
sided power spectral density, 2$^. , with ^n given in (3).* 

With identical receivers in the sum and difference channels which 
amplify the signals (2) by a factor Ko and reduce the center frequency 
from oift to an intermediate frequency, wc = 2Tr X 60 mc, the input to 
the sum and difference channels of the coherent angle-error detector 
can be represented by 

e,,(0 = IU[E,{R,^) cos (o>.-f + did)) 

H- X,{i) cos {(^ii + eiii)) 

-f F,(i)sin (a).-i+ Bid))] 
Cd,(t) = Ko[v^E,(R,0) cos (wit ■^v-\- ei(t)) 

+ Xj(t) cosii^it + d.(t)) 

+ Y^(t) sin (o^it -f- eXt))] 

where X,{t), Xd{t)j Y,{t), Yd(t) have identical one-sided power density 
spectra 2*^ band-limited by the IF bandwidth, Bip , hence have mean- 
square expected values, X- = Y- = ^nBjf . 

III. LINEAR ANALYSIS OF COHERENT ANGLE-ERROR DETECTION SYSTEM 

A block diagram of the coherent angle-error detection system is 
given in Fig. 2. The coherence between the input error signal and the 

* The fluctuation of Si due to Doppler and random phase effects is assumed to 
have negligible effect on the power spectral densities of A'(() and Yit). 
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Fig. 2 — Coherent anglp-enor detection system block diiigram. 

local reference signal is achieved through the action of the phase-lock 
loop. The mixers and detectors indicated by circles in I-'ig. 2 are assumed 
to be ideal multipliers with uuity gain. The ACiC actiou is assumed to 
respond perfectly to variations in the sum channel signal level, so the 
gain of the 60-mc IF amplifiers in both channels can lie expressed as 



K^(e.) = -J-, 



E 



EAH.&)' 



E = constant. 



(6) 



Using (5) and (G), the input to the mixers in each channel in Fig. 2 
can be \\Titten 



e.,(0 - KaE Fcos (o=,t + e.) + Y ^^^ ^"'^ + ^'^ 
+ ^sin(a),f -\-ei)\ 



(7a) 



KoE 7 



edJO = KoE I 1J0 cos {uiit + V + flO + ^ cos (a.,i + 5,-) 
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E, 

(7b) 

■j- ^ sin (oyd -\- Bi) \ 

where the dependence of E, on B and /3 and the time dependence of 
di , X, and Y is understood, but not indicated explicitly in (7) and the 
subsequent analysis, to simplify the equations. 

The other input to the mixers is the output of the voltage-controlled 
oscillator (VCO) in the phase-lock loop, which has the form 

e„(0 = E„ cos (o).,/ + d„) (8) 

where aj„ = 27r X fi5 me and fl„ = d^{i) is the instantaneous phase of 
the VCO output, determined by the operation of the feedback loop in 
the sum channel, which is discussed in the next section. Multiplying 
(7) and (8) and taking only the low-frequency (5-mc) components 
gives at the outputs of the 5-mc IF amplifier in the sum and difference 
channels (see Fig. 2) 

e,,{t) = e\ cos (aj„( - 0i + e„) -1- ^ cos {wot - ei-\- 6,.) 

- ^ sm (coot - e, + 9.) 
ed,{t) = e\ v0 cos (co„( ~ <p ~ Bi -\- e,) 

E, -"" '"""' "' ' ""^ E, 



(9) 



+ ■j=^ cos (c^ot - Bi -I- e„) ~ ~ sin {uoi - Bi+ e,.} 



where 



Es = ^ KdUKmEE,, iC™ = mixer gain, (volts)-! 

Wo = aj„ — 6J,- = 27r X 5 mc. 

The difference channel voltage, Cd^it), is applied to the coherent 
detectors. The other input to these detectors comes from the 5-mc 
reference oscillator whicli produces signals 

fr^ = Er cos a)o( 

(10) 

Cr^ = Ef sin woi 

for the detection of the desired x and y error components given in (1). 
The phase of these signals is the reference phase, 6, = 0. 
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The low-pass filters following the coherent detectors pass only the 
baseband components of the products (eja-e,^) and (eda-er^). Using 
(91 and (10) these baseband components are 



eji 



= a\ ij^ cos i<fi + Oi - B,) + -^^cos (di - $,) 



+ |^^in(0, -fijj 



ej„ = .4 r,8 sin dp -{- d, - 6.) + ^^sin (di - 9,) 



(11) 



where 



A = ^KdEzEr = channel ampHfication factor 



Kd = detector gain, (volts)"^ 

If the phase-lock loop is tracking properly, the phase of the VCO 
output, 6v , will follow closely the phase of the input, Bi . Assuming 
that the rms value of {di — B„) is small compared to 1 radian, then the 
following approximations liold with high probability* 

sin (6i - 6r)= (di - B„) « 1 

cos (Bi - e,) = 1 

and the coherent detector outputs (11) can be expressed in the approxi- 
mate foi-m 



Cdx - Atj^ cos ip - A\v0 sin ip - ^j(^i - ^v) + ^ ^ . 
Cdg =^ AvP sin ip + A(j}fi cos <p + tt )(fi. - 6„) - A 



E, 



(12) 



The first term in each of the expressions in (12) is the desired error 
component, given in (1), amplified by the total difference channel 
gain, At}. The second term represents the perturbation due to the lack 
of perfect phase coherence, while the third term represents the contri- 
bution of thermal noise in the net noise bandwidth of the difference 
channel. 



* The validity of these assumptions is discussed in the next section. 
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The achievement of good phase coherence in the detector outputs (12) 
over the expected range of satellite tracking conditions is the objective 
of the phase-lock loop analysis and design described in the following 
sections. 

IV. PHASE-LOCK LOOP ANALYSIS 

The coherent detection of the control error signals depends on the 
performance of the phase-lock loop in the sum channel (see Fig. 2). 
The loop must be capable of following the change of phase of the input 
signal due to frequency instability of the source and Doppler shift and 
also discrimmate against random phase fluctuations caused by thermal 
noise. These requirements are somewhat contradictory, the former 
requiring a wide loop bandwidth and the latter requiring a narrow 
loop bandwidth. Proper design of the phase-lock loop must therefore 
be based on the best compromise of those requii-ements consistent with 
the expected variation of the signal phase and the expected random 
phase fluctuation. 

The sum channel voltage, Ca/O at the input to the bandpass limiter 
is given in (9). The effect of the limiter can be closely approximated as 
multiplying this voltage by a limiter suppression factor, a, which 
increases from to 1 as the signal-to-noise ratio at the hmiter input 
increases from to « . Tliis hmiter action is discussed in Appendix A. 

The hmiter output voltage, e,^(t) = ae„Jt) is applied to the phase 
detector in the sum channel, together with the reference signal, Cr^ , 
given in (10). The baseband component of the phase detector output 
is therefore 

€b{t) = (er„-Q:e63)bnfleban.l 

or, from (9) and (10) 



eb(t) 



= aA^l^l + sin (e, -e.) -^^ cos (0, - ^,.)] . (13) 

To develop an approximate linear model for the phase-lock loop, the 
following two assumptions are made: 

(i) the phase error (6i — d-^) is sufficiently small to permit the approxi- 
mations sin {di — 9t,) = di — 6„ , cos (0i — do) — I, and 

(ii) the noise component, X(t), is assumed small in the rms sense com- 
pared to the signal amplitude, Es . 

Using these assumptions, the phase detector output (13) can be written 
in the approximate form 

c,U) ^ aAh(t) - 8Al) - ~fj. (14) 
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The low-pass filter passes the baseband component eb{t), producing 
a voltage ef(t) which eaiiscs the frequency of the VCO to vary from 
the center frequency, w, = 2t X fio mc, with a proportionahty constant, 
Kr radians per second per volt. The instantaneous phase of the VCO 
output is therefore 



9v{t) = KAesit) c7/, radians. 



The transfer function of the VCO over the range in which this pro- 
portionality holds is then 



e/(s) 



s 



(15) 



For a phase-lock loop which is stable over a large range of loop gain 
variations and which has zero steady-state phase error for a phase 
ramp input (step frequency cliange) the low-pass filter should have a 
transfer function of the form 



F(s) = 



p/(s) ^ 1 + ts 
Cbis) s 



(16) 



This transfer function has lieen shown^ to yield optimum loop per- 
formance for phase ramp inputs in the presence of white noise, where 
the performance measure is the mean-square error caused by noise 
plus the integrated-squared transient error to tlie ramp input. 

Tlie transfer function (10) can be closely approximated by the 
operational amplifier circuit shown in Fig. 3, which has the transfer 
function 



FM = -, 



1 + R-iCs 



(17) 



1 + iiRiCs 
where — m is the amplifier gain under load without feedback. The 






Zj. » R, 



■AAV 



-u 



U»1 



Fig. 3 — Operational amplifier low-pass filter circuit. 
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assumptions made in deriving (17) are that the input impedance, Z, , 
of the amplifier is very large compared to ffi and that ^ ^ 1 (typically 
10*) in the low-frequency range. With RiC = 1 second and /i of the order 
of 10*, the transfer function (17) reduces to 



Fis) ^ ~ 



1 + R2CS 



(18) 



which is the desired form (16), with R2C = t. The negative sign in 
(18) is incidental provided the sign of the total gain around the loop 
gives negative feedback. 

Using (14), (15), and (IG), the linear equivalent block diagram for 
the phase-lock loop is shown in Fig. 4. The objective of the design is to 
minimize the mean-square value of the random phase error, Bt{i) = 
di{t) — 6„{t), consistent with the requirements on the dynamic tracking 
capability of the loop. 
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Fig. 4 — Block diagram of phase-lock loop based on linear analysis. 
The total transfer function around the loop is 



e.is) 



s~ 



(19) 



where K = AK,. . Ilf R^C is not unity, as was assumed in (18), then 
K = AKr,/R,C.] 

For the analysis of the phase error, 8^ , due to the noise sources 
8n(t) and Ys(t), we let the signal phase 8„ he zero and obtain from 
Fig. 4 the transfer function for 8, in terms of G(s) in (19) 
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Since the one-sided power density spectrum of d„(t) is *» , given in 
(4), and tliat of Ye(t)/E, is 2^nfE^ where ^^r is given in (3), and since 
thpy are uncorrelated random variables, the one-sided power density 
spectrum of the random phase error, de , is 



1 



1 + GO"^) 



2 



1 + G(ja,) 



^^ (21) 



The mean-s(iuare value of this random phase error is then, from 
(19), (21) and integral tables' 



2 



^-l>'-t'^i{^)-U'^-£)-^'-^^^^ 



This expression for the mean-square phase error can be written in 
terms of the undamped natural frequency, ai„ , and the damping ratio, 

f , of the phase-lock loop, as 



4fa)„ r„ E,~ 



(I + 4f^ 



<^/ = T^^— +^-„{^-L_l^) (2.3) 



since «„ = -y/aK, and 2f = ts/oK. 

The proper operation of the phase-lock loop depends upon the magni- 
tude of the phase error remaining less than ■k/2 radians. For the phase 
error due to random fluctuations we can reciuire only that the proba- 
bility of its magnitude exceeding t/2 radians be very small. A criterion 
for this which has been chosen'*'^ as a realistic measure of the threshold 
of the phase-lock loop is that the mean-square value of the total phase 
error be restricted by 

0-/ g i rad^ (24) 

For a normally distributed random phase error with zero mean and 
variance a^ , this criterion implies that the prolmbility of exceeding w/2 
is exceedingly small (about 1()~'). This criterion also gives validity to 
the first assumption made above in obtaining the linear model of the 
phase-lock loop, namely that sin 9, = 6e and cos 9^ — 1. The error in these 
approximations is quite small jirovidod 

\ 6,\ ^ 0.57 radian 

which holds with approximately 90 per cent probability when the 
(condition (24) is satisfied. 

The second assumption made above for the linear model is not strictly 
justified in the region of threshold, where the signal-to-noise ratio in the 
3-kc bandwidth at the phase detector input will typically be less than 
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unity. However, comparison studies with a digital computer simulation" 
which was implemented for the angle-error detector in the precision 
tracker system have indicated that the linear model estimate of the 
mean-square phase error (23) is sufficiently accurate even in the vicinity 
of threshold to justify its use for the analytical design optimization of 
the phase-lock loop. The digital computer baseband model of the phase- 
lock loop includes the in-phase noise term as well as the quadrature noise 
term in (13), and the sine and cosine operations of the phase detector 
(see Ref. 0). A comparison of the computer simulation data with the 
linear analysis data for the design examples considered in Section VI is 
given at the end of that section. 

An important parameter in the phase-lock loop analysis is the effective 
noise bandwidth, B^ , of the loop, defined i>y 



B, 



=/: 



df. 



Br. = 



'o 1 1 + G(s) 

The loop noise bandwidth for the system under consideration has 
already been evaluated in the second terms of (22) and (23), namely 

which increases with the hmiter suppression factor a and hence in- 
creases with the signal -to -noise power ratio at the limiter input. This is 
the desired adaptive feature which the limiter provides in the phase-lock 
loop operation, since for small S/N (long-range condition) the loop 
bandwidth is small, decreasing the mean-square error due to thermal 
noise, while for large S/N (short-range condition) the loop band- 
width is large, providing improved phase tracking accuracy for the 
greater Doppler frequency rate of change occurring at short range. 

The Doppler frequency variation as a function of time is approxi- 
mated by a frequency "ramp" input having a constant slope of magni- 
tude U3 for a duration T,i , as indicated in Mg. f). Since this approximate 
function represents a somewhat more difficult variation for the loop to 
track than the actual Doppler variation, tlie evaluation of the loop 
tracking accuracy and transient behavior based on the approximate 
input should serve as a consei'vative basis for design. 

Neglecting the thermal noise terms in (13), the phase detector 
output reduces to 

eb{t) = aA sin [di(t) - Ojt)] volts (26) 

where d{(t) = ^cht^ is the phase input corresponding to the frequency 
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Fig. 5 — Typiciil Doppler freqiieiirv- vnritilion and pircowisc linear approxiination. 

"ramp" input discussed abovo. Since this output is l)Ounded in magni- 
tude Ijy ±A volts, the low-pass filter which follows the phase detector 
should have a linear input dynamic range of at least this magnitude. 
In addition, the voltage-controlled oscillator should have a linear 
frecjuency range at least as large as the expected maximum Doppler 
shift. With these two design requirements satisfied, the only essentially 
nonlinear element in the phase-lock loop circuit is the phase detector 
(see ]''ig. 6a). 

An analysis of the response of this nonlinear circuit to a frequency 
"ramp" input is given in Appendix B. h'rom this analysis it is concluded 
that for adequate phase-lock tracking of the Doppler shift, the loop 
gain should satisfy the condition 

aK > 2tin.a. (27) 

where a),n.ix is the maximum Doppler rate in rad/sec". 

When condition (27) is satisfied, the steady-state phase error due to 
tins Doppler rate will not exceed tt/Q radian (sec Appendix B), and the 
loop i-esponse will closely approximate that of the linear second-order 
circuit, shown in Fig. 6(h). Por this circuit with the input 8i(t) = ^ut^, 
the Laplace transform of the phase error is 



eM = 



s[s~ + 2fa)„s -I- oi„^] 



where, as before, w„ = s/aK, 2f = r^/aK, and K = AK^ . The time 
response of this type of liTiear second-order system is thoroughly dis- 
cussed in elementary texts on linear circuits or control system theory.* 



* Sec, for example, Ref. 10. 
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Fig. 6 — Phase-lock loop circuits for Doppler-shift analysis: (a) equivalent 
circuit including nonlinearity of phase detector; (li) approximate linear circuit, 
valid for | ffg | S tt/G radian. 

The two per cent settling time* for the transient of the phase error is 

approximately 



8f 
Ta = — = 4t seconds 

0}„ 



(28) 



for the overdamped case (f > 1). This also serves as a good approxi- 
mation of the duration of the transient for the underdamped case in 
the range 0.7 < f < 1. 



V. PHASE-LOCK LOOP DESIGN 



The design procedure which will be followed is first to ascertain the 
system requirements necessary for minimum acceptable performance 
and second to optimize the performance within the range of parameter 
adjustment available to the designer. 



* Defined as the time required for the transient response to settle down to 
within two per cent of the steady -state value. See Ref. 10. 
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The primary factors external to the phase-lock loop which affect its 
performance capability are : 

(i) The equivalent coherence time, tc, , of the input signal, which 
characterizes the random frequency fluctuations of the various system 
oscillators which affect the instantaneous frequency of the signal into 
the phase-lock loop. The coherence time should be made as large as 
possible, but is primarily limited by the frequency stabihty achievable 
in the small satellite transmitter and hence will be considered a fixed 
parameter not available for phase-lock loop design adjustment. 

(u) The noise-to-signal power ratio at the input to the phase-lock 
loop. Since this ratio varies with the range of the satellite and with the 
IF bandwidth preceding the loop, it is desirable to characterize the 
relative "noisiness" of the system independent of range and bandwidth 
variations. The thermal noise power is given by 

AT = $j^B,F = kT.gBiP watts 

while the average signal power is^'-'^ 

where 

Pt = transmitter power, watts, 
Gt = transmitter antenna gain, 
R = transmitter to receiver range, and 
Ar = effective area of receiving antenna (same units as R"^). 

The noise-to-signal power ratio can then be expressed as 

The factors in parenthesis in (29) are constant* characteristics of 
the satellite-ground system and will be represented by a single constant 
called the "receiver noise index," denoted by the symbol K , and 
having the units of seconds/ (distance) . Then, 

S E. 
The noise index, ^■, , is taken as the second fixed parameter! of the 

* The transmitter antenna gain, Gt , will not actually be constant unless the 
radiation pattern is uniform or tlie satellite is properly attitude controlled. It is 
assumed essentially constant in this study. 

t From (30) it is apparent that the receiver noise index, k, , corresponds to 
the unit bandwidth noise-to-signal power ratio in the receiver when the satellite 
is at unit distance from the receiver. 



-BiP = hR'B,^. (30) 
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external system affecting loop design. 

(in) The satellite orbit characteristics.* The orbit characteristics 
affect the performance capability of the phase-look loop as follows: 
First, the maximum and minimum communication range will determine 
the variation of tlie noise-to-signal ratio as shown in (30). Second, the 
maximum range rate, Rum^ , will determine the maximum Doppler 
shift 

Ao) I ,„»^ = — — rad/sec (31) 

c 

where 

uib = 2t X satellite beacon frerinency, and 
c = velocity of light. 

Finally, the maximum range acceleration, R,,,^^ , will determine the 
maximum Doppler rate 

tCft/tniflK 1 / 2 /nn\ 

o},n:iK = rad/sec . (32; 

c 

Since the maximum w occurs at the range of closest approach for all 
possible satellite passes, the limiter suppression factor a will be at its 
maximum value, giving the largest loop gain, aj„B.xK- From (27), the 
basic lower limit on the fixed loop gain constant, K, is then 

h > sec . (33) 

The total mean-stjuare error due to random fluctuations in the system 
is given by (22) or (23). Using (25) and (30), the mean-square error 
can be expressed as 

<r' = ^-^^^^ + KR'B, , VBd' . (34) 

The second term in (34), due to thermal noise in the system charac- 
teiized by /iv , increases with the square of the range. This increase is 
somewhat offset by the reduction of the loop bandwidth, Bl , which 
decreases approximately as the first power of the range. This reduction 
of Bl , howeve]-, increases the mean-square error due to random pha.se 
fluctuations of the signal, given by the first term in (34). Hence, the 
total mean-square erroi" will be maximum at the longest range condition. 

* Derivation of pertinent orbit characteristics is given in Appendix C. 
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From (24), a requirement which the loop design must satisfy is there- 
fore given by 

(ff/)ff=/r,„„^ ^ 8 radl (35) 

As a function of Bi, , the minimum possible value of a^^ is, from (34) 



mm ai = 
when 



''[ti'^w)] '''' 



-'^"p-'4(^^^y- 



Tt is apparent from (36) that unless 



J-rR^ < i- (:J8) 



Tee 



32 



the (rendition (35) cannot be satisfied for any values of loop t)andwidt.li 
and damping ratio, or cquivalently for any values of loop gain, aK, 
and time constant, r. This corresponds to the minimum signal power 
condition given by Develet. 

The two basic design requirements which the external system imposes 
upon tlie phase-lock loop design arc given by the inequalities (33) and 
(38). The first requirement does not appear critical, since for satellite 
communication systems the value of Wnmx is unlikely to exceed about 
2 X 10^ rad/sec^ (see Appendix C). Since a^ax - 1, condition (33) 
requires that 

/v > 4 X 10* sec^" 

while values of K on the order of ten times this lower limit are achievable 
in present phase-lock loop designs. 

The second requirement, (38), is more critical, since it depends 
entirely upon the hxed parameters of the external system. If the in- 
e(|uality (38) is not satisfied, either the system noise index must be 
decreased or the coherence time of the satellite transmitter must be 
increased before satisfactory operation of the phase-lock loop can be 
achieved at maximum range. When this letjuirement is satisfied, the 
optimum loop design with mspect of the total mean-square error is 
achieved by making the loop bandwidth equal to the optimum value, 
given by (37), and by making the damping ratio, f, as large as is con- 
sistent with satisfactory transient response of the loop. 

Fig. 7 depicts in graphical form the design requirements discussed 
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above. For a given value of 



Tee, = 



Toe 



1 + (1/4^) 



(39) 



the contour of a^' = \ rad" is given as a function of KR^ and the loop 
bandwidth Bl ■ Within this contour the mean-square error is less than 
the threshold value of \ rad' and has its minimum possible value for a 
given ktR' on the dashed line associated with each contour. The vertical 
line defined by krR' = KR^ax must intersect the contour for the ap- 
propriate value of Tee in order for satisfactory loop design to be achieved. 
If this critical requirement is satisfied, then the optimum design is 
given by adjusting the loop bandwidth, Bl , to equal the value obtained 
from the dashed line within the contour at the particular value of 

krR\ 

The input-adaptive adjustment of B^ by means of the limiter sup- 



4000 ■ 




,0-6 <: " •> =10-5 ^ " " "lo-" ^ * » " 10-3 = 

RECEIVER UNIT BANDWIDTH NOISE-TO-SIGNAL POWER RATIO, krfif IN (CPSj' 



Fig. 7 — Phase-lock loop threshold contoiirB. Inside the contour for a par- 
ticular Tc' the mean-Bquiire phiiBC error, ff,^, is less than the threshold value of 
0.125 rad*, and has its minimum value for a given krR' on the dashed line as- 
sociated with the contour. Indicated on the ahscissa are the expected maximum 
values of krR^ for the autotrack and the precision tracker. 
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pression factor, a, can approximate this optimum adjustment of Bl . 
The optimum Bl , given by (37), varies inversely with the range, R. 
The actual value of Bl , given by (25), varies directly with a, but a 
itself varies inversely with R, as shown in Appendix A. 

Before the optimum bandwidth can be determined from (37) or 
Fig. 7, the value of f at maximum range must be chosen. As was pointed 
out at the end of Section IV, the settling time of the loop is approxi- 
mately 4t seconds. If Tm denotes the maximum tolerable settling time, 
then this requires that 4r ^ Tm . This also places an upper hmit on f 
when the loop bandwidth, Bl , is fixed at the optimum value (37), 
since, from (25), Bl and f are related to t as follows 

Bl=^(H"--{-1). 
4t 

To satisfy the maximum settling time restriction, then 

4r = ^^-^ ^ Tm (40) 

Bl 

which imphes that f cannot be arbitrarily increased while holding Bl 
fixed at the optimum Bl . Using (37), (39) and (40), the upper bound 
on f when Bl = opt Bl can be expressed in terms of Tm and the external 
system parameters as 



^2 < ?■ 2 = 1 



T„ + 1 p _ 1 



2t,XR' 



(41) 



The lower bound on f when Bl = opt Bl follows from the require- 
ment that niin it/ < \ rad''. Using (3fi) and some algebraic manipulation 

f 2 > r ' = ^^"^' (42) 

which is a finite real lower bound on f only when the basic requirement 
(38) is satisfied. 

Since f varies with range, the upper and lower bounds given above 
should be evaluated at the same range, R. It is apparent from (41) 
and (42) that the least upper bound and the greatest lower bound on f 
both occur ni R = Rn,B.x , and that the bounds constrict the range of f 
as the noise index, l\ , increases. As fcr decreases, the limits separate, 
allowing a wide range of f . Howe^-er, from the point of view of relative 
stability and fast transient response, f should not be less than 0.7 nor 
greater than about two ; also, from the point of view of minimizing o-f , 
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it can be seen from (1^6) that there is negUgible improvement in in- 
creasing f beyond about two. These restrictions on the range of f in 
the optimum design of the phase-lock loop can be summarized by 

max (0.7, U) < ^ ^ min {2.0, U) (43) 

where ^n and f™ are the upper and lower bounds defined in (41) and 
(42), respectively. 

Assuming that the set of positive real values of f satisfying (43) is 
not empty, the optimum choice of f (with respect to minimizing the 
mean-square phase error) is the least upper bound value given in (43). 
Using this value, the optimum value of the loop noise bandwidth, Bi, , 
at maximum range is determined from (37) or Fig. 7. 

The selection of the optimum f and Bl at maximum range fixes the 
value of the loop filter time constant, r = R2C , as 

^ 4ropt' + 1 ^ /7'„/4 , .c.„t - r.» 

^ 4(5.)o„t 4.25 (44) 

and the value of loop gain at maximum range as: 

a,ni„ K = %^ . (45) 

When the IF bandwidth is specified, a,nin is determined from Fig. 
11 in Appendix A, with R = i?,nnx • lOiowingamin , the loop gain constant, 
K, is then determined from (45). Since K has the lower bound given by 
(33) and certainly an upper bound dictated by practical equipment 
considerations, the range of a may have to be controlled through the 
selection of Bu- to give a value of K which is compatible with these 
bounds. 

VI. TELSTAR SYSTEM DESIGN EXAMPLES 

To illustrate this design approach, sample designs will be considered 
for both the precision tracker and the autotrack systems for the Telstar 
experimental program. Based on analysis of the expected orbit (see Ap- 
pendix C), and on preliminary system data, the system parameters as- 
sumed for the design examples are given in Table I. While the IF band- 
width values are not necessarily fixed, the 200-kc bandwidth given in 
Table I for the 5-mc channels (see Fig. 2) is desirable from practical 
considerations. 

Considering first the dynamic tracking capability in the absence of 
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Table I — Assumed System Parameters for Design Examples 

Equivalent coherence time, Tc, = 0.02 sec 
Receivei' noise index (see Section V) 

PT system Av = 1.5 X 10^^ sec/knm (kilo-nautical miles) 

AT system Av = 2.0 X lO"* sec/knm 

Orbital data: 

Maximum range, i?mni = 5 knm 

Minimum range, Rmn, == 0.5 knm 

Maximum Doppler shift, | V I — 100 kc 

Maximum rate of .shift, | w | = 5620 rad/sec*, at 0.5 knm 
IF bimdwidth preceding litniter, B/p = 200 kc (both systems) 

Nominal range of loop gain constant, K = 10^ to 10' sec~= 

Maximum tolerable settling time, Tm = 0.1 sec 

noise, the dynamic range of the voltage-controlled oscillator in the 
phase-lock loop should lie about ±150 kc to avoid saturation effects 
when the maximum Doppler shift is ±100 kc, and to allow for some 
drift of the center frequency during operation. 

The I'ange of variation of the limiter suppression factor, a, can be 
determined from Fig. 11 in Appendix A, using the values of k^ , R, 
and BjF in Table I 

Precision Auto- 

Tracker track 

Max. range, a = a,nh,- 0.10 0.28 

Min. range, a = a,^^^' 0.79 0.97 

Since the maximum Doppler rate occurs at minimum range, the con- 
dition (33) that the maxinuun steady-state error be less than x/6 
radian requires that tlie loop gain constant, K, satisfy 

Precision 
Traclier Autotraclt 

K > '^ ^ ''"'^^*' = 1.42 X 10' 1.15 X 10^ 

Both of these lower limits ai'C well below the lower nominal value of 
10^ hoc"' given in Table I. Using this lower nominal value the maxi- 
mum steady-state error in tracking the Doppler shift will be less than 
0.08 radian for both systems [see (50), Appendix B]. 

The upper and lower bounds on f at maximum range condition, using 
(43) and Table I, are 



Precision 
Tracker 



o.fiil fi.7r, 

0.70 ^ ^ ^ 12.0 



Autottack 



0.15\ ^ ^ ^ [2.95 



0.70J ^ = \2.0 
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Taking the least upper bound as the optimum value gives aiR = Rmax 

Precision I 

Tracker Aulotrack 

r^pt =1.76 I fopt = 2.0 

The min ct/ and optimum Bi, at maximum range, using (36) and (37) 
and the above values of fopt , are 



min (T^ = 0.1 rad 
o-ptBL= 135 cps 



min ct/ - 0.036 rad^ 
opt Sl= 364 cps. 



Finally, using (44) and (45) and the values of amin given above for the 
200-kc IF bandwidth, the optimum values for the phase-lock loop 
constants t and K are 



T = 0.025 sec 
K = 2 X 10' s 



T = 0.012 sec 

ff = 4 X 10' sec"'. 



Since both values of K are greater than the lower bounds given above 
and are within the nominal range given in Table I, the desired 200-kc 
bandwidth need not be changed. 

The performance of the systems as a function of the range, R, using 
thet-e design parameters, is .summarized by the curves shown in Figs. 8 
and 9. These curves show how the input adaptive loop gain, aK, helps 
provide near-optimum design at ranges other than the maximum range 
where the optimum design was accomplished. 

The mean-square phase error obtained from this linear analysis was 
compared with values obtained from the more accurate digital computer 
simulation of the phase-lock loop described in a separate paper." Results 
for the two design examples considered above, which are shown in Fig. 
10, demonstrate that although the error in the linear model estimate in- 
creased somewhat as the signal-to-noise ratio decreases, there is no 
drastic breakdown in the accuracy of the linear model. For these and 
two other designs tested, the linear analysis estimate of the threshold 
signal-to-noise ratio was within 1.5 db of the digital computer results. 

VII. CONCLUDING REMARKS 

The material presented in this paper was part of a design study 
conducted during the initial construction of the ground station tracking 
systems for the Telstar program. At that time there was some concern 
about the capability of the angle-error detector to maintain phase-lock 
at the longer satellite ranges, due to the small beacon signal power 
available and the uncertainty about the coherence time of the 4080-mc 
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Fig. 8 — Performance curves for precision tracker (PT) optimum design 
example: (a) phase-lock loop parameter variations with range from ground 
station to satellite; (b) phase-lock loop performance measures as a function of 
range. 

beacon signal. This was particularly critical in the detection system for 
the precision tracker because of the significantly lower gain of the 
precision tracker antenna compared to the horn-reflector antenna. For 
this reason, the means for achieving optimum design of the phase-lock 
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Fig. 9 — Performance curves for auto track optimum design example: (a) phiise- 
lock loop parameter variations with range from ground station to satellite; (b) 
phase-lock loop performance measures as a function of range. 

loop at maximum range was a crucial consideration in the initial design. 
One of the fortunate results contributing to the highly successful 
operation of the first Telstar satellite experiment was the excellent 
phase stability achieved in the satellite beacon transmitter. Measure- 
ments of the mean-square phase error under strong-signal conditions 
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Fig. 10 — Compariaon of metin -square error predicted by linear analysis with 
results from digital computer simulation of phase-lock loop, using design values 
from the two examples in Section VI. 



(where the major contribution to the error is due to phase jitter in the 
beacon signal) indicated an effective coherence time of about 0.1 second 
instead of the 0.02-second value assumed for the design examples 
above. This higher coherence time makes the noise bandwidth and 
dampmg ratio adjustment in the phase-lock loop much less critical, as 
can be seen from the threshold contours in Fig. 7. 



1966 THE BELL SYSTEM TECHNICAL JOURNAL, SEPTEMBER 1963 

More detailed descriptions of the final design and the performance of 
the tracking systems in the first Telstar experiments are given in a 
series of papers ^•'^•^ appearing in an earfier issue of this journal. 

APPENDIX A 

Effective Gain of Bandpass Limiter 

Two factors which characterize the operation of an ideal bandpass 
limiter are: 

(i) The total power output of the limiter remains constant/^ i.e. 

.So -^ No = C (46) 

where So = output signal power, A^o = output noise power. 

(ii) When a sinusoidal signal and narrow-band Gaussian noise are 
apphed to the input of a bandpass limiter, the output signal to noise 
ratio is related to the input signal to noise ratio by^^ 

^ = X ^ (47) 

where the factor X, given in Fig. 5 of Ref. 13, varies from t/4 to 2 as 
the input signal to noise ratio varies from zero to infinity. 

When no noise is present, we assume that the output signal power 
equals the input signal power (any fixed gain in the bandpass limiter is 
absorbed into the loop gain constant, K). Then, from (46), when 
A^o - 

Sq = C = Si 

so that, when noise is present 

So+N,= Si . 

A little algebraic manipulation of this expression gives 

So 

-So A^o 



Si 1 I ^ 



Using (47), we obtain 
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So 

Si 






1 +x 



Si 
Ni 



(48) 



where a is called the limiter suppression factor." 

Thus, the limiter has the effect of reducing the signal power from S 
to aS, and as a consequence the effective loop gain is reduced from K 
to aK. The factor a defined in (48) varies from to 1 as the input 
signal-to-noise ratio varies from to x. From (30) in Section V, the 
input signal-to-noise ratio can be expressed as 

Using this expression in (48) gives 



[-^^r' 



« = 1 1 + 

which shows the inverse dependence of a on the sateUite range, R. 
With values of \ obtained from Fig. 5 of Ref. 13, the hmiter suppression 
factor, a, is plotted as a function of the input signal-to-noise ratio in 
Fig. 11. 




LIMITER INPUT SIGNAL-TO-NOISE RATIO, rr = (KpB^Bip) ■ 



rig. 11 ■ — Limiter suppression factor, a, as a function limiter input signal-to- 
noise ratio. 
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APPENDIX B 

Nonlinear Analysis of the Phase-Lock Loop 

The phase-lock loop equivalent circuit, which includes the sine- 
function nonhnearity of the phase detector, is shown in Fig. 6(a). In 
terms of the phase error, 6e , and its time derivative, We , the differential 
equations governing the response of this circuit to the frequency-ramp 

input, toi = 6it, are 

dde 

(49) 



dr = ^" 



-j-^ = w — 2fwnaj„ cos 6e — w„ sin dg 

(JiL 

where, as in Section IV, we define u„ = aK, 2f ^ rwn . 

The values of Qe and w^ which satisfy the equilibrium condition that 
the right-hand side of (49) vanish are 



(e.)e^ = sin-^ {^\ , UU = 0. (50) 

If the phase-lock loop "locks-on" to the frequency ramp input, the 
frequency error is zero, but there is a steady-state phase error given by 
equilibrium value in (50). A necessary condition for the existence of 
this phase-locked response is that 

aK > ci) 

i.e., the total loop gain must exceed the input frequency rate. 

To analyze further the response of this circuit, it is convenient to 
normalize (49) in time and frequency with respect to the parameter 
a)„ . Defining the symbols 

V = U>e/oin 

T = ti)/a)„ = (jsfaK. 
the differential equation (49) can be written. 

de ^ 

dx 

(51) 

-— = r — 2fv cos d ~ sm9 
ax 

where 6 = Beit) is assumed in the remainder of Appendix B. 
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The solutions of these equations for given initial conditions and 
given value of normalized input rate, r, describe trajectories in the 
normalized (6, v) state space. The slope of these trajectories in this 
state space is, from (51) 



dp 
dd 



r — sm 



- 2^ cos e. 



(52) 



For initial conditions ^ = 0. j' = (corresponding to the circuit 
being in steady-state phase lock with constant frequency input prior 
to the onset of the frequency-ramp input), two sets of trajectories 
obtained from numerical integration of (52) are shown in Fig. 12. In 
Fig. 12(a) a value of r = 0.966 is chosen to illustrate the response when 
the loop gain exceeds the frequency rate only slightly. It can be seen 
that for f ^ 1, the phase error tends to the steady-state value of 1.31 
radians (= sin"' 0.960) with only small overshoot. For f = 0.707, 
however, there is a large overshoot which actually exceeds jr/2 radians, 
but eventually i-oturns to the steady-state value; for f < 0.707 the 
phase error docs not reach the steady-state: i.e., the circuit is unable to 
"lock-on." 
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Fig, 12 — Trnjectories of the phase-lock loop response to frequency ramp 
input for viirioiis viiUies of the damping ratio; (a) input frequency rate, i, nearly 
eiiniil lo tlio loop gain, aK; (b) input frequency rate, (I*, equal to one-half the 
loop gain, aK. 
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To avoid large steady-state phase errors and large peak phase errors 
with the attendant likelihood of random perturbations causing the 
circuit to fall out of phase lock, the following conditions should be 
imposed 



(53) 



Fig. 12(b) shows a set of trajectories with conditions (53) satisfied. 
The response closely approximates that of the linear second-order 
system obtained by letting sin 8=6, cos 9 = 1 in (51). 

APPENDIX c 

Satellite Orbit Characteristics 

The parameters of the satellite orbit which affect the design of the 
phase-lock loop are evaluated in this appendix. The effect of the oblate- 
ness of the earth and other perturbations upon the satellite orbit is 
neglected in this analysis. However, since this effect does cause rotation 
of the perigee, the maximum range, minimum range and maximum 
Doppler effects are derived considering all possible locations of the 
perigee relative to the ground station. 

The geometry and terminology of the analysis is shown in Fig. 13(a). 
It is sufficient to consider only the condition when the ground station 
is in the plane of the orbit in order to derive all the parameters needed. 

(i) Minimum and Maximum Communicating Range. It is obvious that 
the mininunn possible range occurs when the satellite passes overhead 
at perigee. Therefore 

Rnua = Rp = perigee altitude. (54) 

From the point of view of the satellite, the maximum possible range 
to any visible point on earth occurs when the satellite is at apogee and 
the range is taken along the tangent to the earth's surface. The satellite 
would then appear on the horizon at maximum range to a tracking 
station anywhere along the locus of these points of tangency. 

Since, however, the satellite must be at a small angle, tpk , above the 
horizon before communication is feasible at maximum range, the 
conditions at the maximum possible communication range are as shown 
in Fig. 13(b). In terms of the angle <pc in Fig. 13(b) 

flmax^ = r^"* + r^ - 2rarQ cos ip^ 
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where 

To = radius of earth 

Ta = Ra-l- ro, Ra = apogce altitude. 
Now for ifih small (less than about 10°), the angle <pc is very closely given 




PERIGEE 
ALTITUDE 



(a) 



SATELLITE 



TRACKER 




Fig. 13 *— Satellite orbit diagrams; (a) geometry and terminology for satellite 
orbit; (b) conditions for maximum communication range. 
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by 

(Pc — ipi — fh, vi = COS - 1 — 



■fe)- 



Therefore, in terms of known parameters and a given horizon angle 
iph , the maximum possible communication range is 

ffn.»> = Va + r^ - 2r„ro cos {ipt - <Ph)f (55) 

where 

ipt = cos~^ (ro/ra). 

(it) Maximum Doppler Shift and Rate. The orbit parameters needed 
to determine these Doppler effects are the maximum range rate, R,nax , 
and the maximum range acceleration, flmax ■ The parametric equation 
of the satellite orbit in polar form corresponding to the choice of coor- 
dinates in Fig. 13(a) is 



where 



2r„r^ 
n„ = 



Ta + r^ 



€ = 



' (1 ' p 



= eccentricity of orbit 



fa + Tj, 

Ta = Ra-\- To 

Tp = Rj,-\- ro. 
Furthermore, from the "law of areas" for motion in a central force field 

^ ^ - rad/sec (57) 

where 

k = GMrm = gro r„ 
G = universal gravitational constant 
M = mass of eartli 
g = acceleration due to gravity at surface of earth. 

From Fig. 13(a), the range, R, for a tracker at angle tp is related to 
the orbit variables, r and 6, by 
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i?^ = r' + ro' - 2ror cos (6 - <p). (58) 

Differentiating this expression and using (57) gives for the range rate R 

R = L[r -n cos {d-<p[+ ^^ sin {d - <p) (59) 

R iir 

and for the range acceleration, R 

R=llr- r, cos (9 - ^)\ + '^^^ + ^^ cos {6 - <p). (60) 

These expressions depend upon r, r, and f, which are determined as a 
function of 6 by the orbit equations (56) and (57). The evahiation of 
^m.« is rather tedious and is most easily obtained for a given orbit by 
machine or graphical computation. It was evaluated for the expected 
Telstar satellite orliit using a part graphical and part analytical com- 
putation, with the results given at the end of this appendix. 

Tiie evaluation of Rmax is quite easily obtained, however, since it 
occurs for the conditions 6 ^ ,p = 0; i.e., when the sateUite passes over- 
head at perigee. The maximum range acceleration is given by 

which occurs when R = RmiK = Rp ■ 

The maxinuun Dopplcr rate varies directly with the maximum range 
acceleration 

.„.,=?^'«,„.^2.A^[^^.(|-l)] (62) 

uhi'i'c 

Ji = satellite beacon freqviency 

c ^ velocity of hght. 

To estimate the luaximum Dopplev rate which might be expected for 
practical communication satelUte systems, we take the following 
conditions as representing practical extremes from the point of view of 
good eonimimication and satellite lifetime 

maximum /& = 5 X 10^ cyc/sec 

minimum perigee, Rp ^ 0.2 knm 

maximum apogee, Ra = 5.0 knm. 
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Using these values in (62) gives as an estimate for the maximum 
expected Doppler rate 

max (wniax) f^ 2 X 10* rad/sec^. 

(m) Numerical Values for the Expected Telstar Satellite. The constants 
needed for the range and Doppler calculations are: 

Rp = 0.5 knni (perigee) 

Ra = 3.0 knm (apogee) 

To = 3.44 knm 
r„, ^ 4.88 knm 
g/c = 3.27 X 10"' sec"'; c = 162 knm/sec 

fb = 4.08 X lO^cyc/sec 

tph = 7.5° (acquisition angle above horizon). 

Using these numerical constants in (54), (55), and (62) gives for 
^min , ^max , aud 6}„,ax the vahics : 

Riuh, = 0.5 knm 

ftmax = 5.0 knm 

w^ax = 5.62 X 10' rad/secl 

The maximum value of R in (59) for the Telstar satellite orbit was 
found to occur when 6 = 340°, ^p = 10° and has a magnitude 

I R Uax = 4 X lO"' knm/sec. 

The maximum Doppler shift is then given by 

|A/1,„,, =^-\R\u... = 100 kc. 
c 
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